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ABSTRACT 

Star-forming galaxies are predicted to contribute considerably to the cosmic gamma- 
ray background (CGB) as they are confirmed 7-ray emitters and are the most numer- 
ous population of 7-ray sources, although individually faint. Even though the Fermi 
Gamma-ray Space Telescope will be able to resolve few star-forming galaxies individu- 
ally, their fractional contribution to the CGB should become far more significant than 
it was for past measurements of the CGB as many of the brighter, formerly unre- 
solved sources will be resolved out. Thus, the clustering feature of galaxies imprinted 
on the CGB might be detectable by Fermi. In anticipation of such measurements, we 
calculate the predicted angular auto-power and cross-power spectra of the CGB from 
normal galaxies. We find that the amplitude of the auto-power spectrum is smaller 
than that for other sources such as blazars and dark-matter annihilation; the shape is 
also characteristic. We also show that the cross-power spectrum with galaxy surveys 
features larger amplitude. Fermi should be able to detect the correlation signature in 
both the auto-power and cross-power spectra at angular scales of ~1-10° after 5-yr 
of operation. Such a detection would be valuable in confirming the level of the star- 
forming galaxy contribution to the CGB, and more importantly, in serving as a tool 
in the effort to discriminate between possible origins of the CGB. 

Key words: gamma rays: theory — large-scale structure of Universe — galaxies: 
evolution — cosmology: theory. 



1 INTRODUCTION 

Star-forming galaxies are confirmed 7-ray sources. They 
emit 7-rays produced in hadronic interactions between 
cosmic-ray nuclei and interstellar gas, and in leptonic 
interactions between cosmic-ray electrons an d sec- 
ondaries with interstellar gas a nd li gh t (e.g. ISteckeil 



197d. Il97j: iFichtel fc KniffenI Il984l; iDermed [iQsl 



Strong. Moskalenko. fc Rcimor 2000l ) . Diffuse 



emission 



from the Milky Way is, in fact, the bright est feature of 
the 7-ray sky, as dem onstrate d by SAS-2 jKniffen et "al] 
1 19731 ). COS-B iMaver- Hassclwa nder et"aLl Il982l ). the En- 
ergetic Gamma-Ray Experimental Telescope (EGRET) 
onboard the Comp ton Gamma-Ray Observatory (CGRO) 
(|Hunter et al.l[l997l ). and by the first-light res ults of Fermi 
Gamma-ray Space Telescope (|Ritz et al.ll2009l ). Other than 
the Milky Way, the only star-forming gal axy detected in 
gamm a-rays is the Large Magellanic Cloud (jHartman et al.l 
I1999I ). because normal star- forming galaxies are individ- 
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ually faint in 7-rays (jPavlidou fc Fieldi l200l[ ). However 
star-forming galaxies are very numerous, and their collec- 
tive emission is likely to make a substantial contribution 
(Pavhdou & Fields 2002, hereafter PF02) to the cosmic 
gamma-ray backgroun d (CGB), measured with EGRET 
jSreekumar et al.ill998h . 

The Large Area Telescope (LAT) onboard Fermi will 
further refine the CGB measurement with improved en- 
ergy and angular resolutions, whereas it will only de- 
tect no more than three additional galaxies as individual 
sources (Small Magellan ic Cloud, M 31, and maybe M 33; 
IPavlidou fc Fields! l200l[ l. Therefore, normal galaxies would 
be a guaranteed source of the CGB for Fermi-hNT , and their 
contribution would be at essentially the same level as it was 
for EGRET. Other contributors such as blazars, on the other 
hand, will be substantially reduced with res pect to their f rac- 
tional contributions to the EGRET CGB l|Steckeilll999h . as 
the Fermi-LAT will resolve many of th em (>1000, depend- 
ing on the luminosit y function; see, e.g.. lNarumoto fc Totanil 
l2006l : lDerm"eill2007l ). Normal galaxies also have a character- 
istic spectral feature — a peak, tracing the hadronic origin 
of their emission. As a result, when the (more spectrally 
featureless) blazar contribution is reduced, the contribution 
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from normal galaxies to the CGB could be dominant at en- 
ergies around a few hundred MeV (PF02)E1 which allows for 
an almost contamination- free set of photons. 

Galaxies are clustered following the large-scale mat- 
ter distribution in the Universe, and this clustering fea- 
ture should be imprinted on the CGB. The anisotropy of 
the CGB has recently been studied theoretically by a num- 
ber of authors, in order to look for signatures of vari ous 
contributing so urces such as blazars (lAndo et al. 2007al lbh ■ 
galaxy clusters ("Ando et al."2007a':'Miniat i et al.ll2007l '). tvpe 
la supernovac (Zhang & Bcacom 2004) , and d ark-matter 



annihilation (lAndo fc Komatsul 120061: lAndo et all l2007bl: 
Cuoco et all | 2007l; iHoop er fc Serpicol 120071: ICuoco et aP 



20081 : ISieeal- Gaskins 200 81: iLee. Ando. fc Kamionkowskil 
2009 : iTaoso et al.i 12008 : iFornasa et al.l l2009l : lAndd l2009l ). 



The same approach should also be taken for the normal star- 
forming galaxies. Should this signature be detected in the 
_Fermi-LAT data, it would be extremely useful in a variety 
of ways: 

(i) As a consistency check. — If an energy range is identi- 
fied spectrally where the normal-galaxy contribution to the 
CGB is believed to be strongly dominant, then the CGB 
photons in this range must exhibit anisotropy properties 
consistent with our understanding of normal-galaxy clus- 
tering. 

(ii) As a powerful tool to disentangle multiple CGB 
components. — Instead of having the CGB strongly domi- 
nated by normal galaxies in some energy range, an equally 
likely scenario is to have a balanced mixture of normal 
galaxy photons and photons from different source classes. 
In this much information as possible is needed to 
disentangle the different CGB contributions. In this con- 
text, the angular power spectrum is as important a clue as 
the shape of the energy spectrum of the contributions from 
differ ent populations (see, e.g., Siegal-Gaskins & Pavlidou' 
I2OO9I on how the two can be combined when information 
from the energy spectrum alone is insufficient to break the 
degeneracy between different components). The importance 
of the angular power spectrum in the case of normal galax- 
ies is further emphasized because their clustering proper- 
ties are very we l l constrained through galaxy surveys ( e . g . , 
ICole et al.|[2005l : iMaller et al.ll2005l : rPercival et al.ll2007l ). 

(iii) As a complement to anisotropy studies of other 
populations. — As normal galaxies provide a guaranteed con- 
tribution to the CGB for Fermi-LAT, the CGB anisotropy 
due to normal galaxies is also a guaranteed background to 
any anisotropics studies using the diffuse background. For 
this reason, it is important to calculate the anisotropy prop- 
erties of the normal galaxy signal and understand its uncer- 



^ Note that while PF02 found that the normal galaxies were likely 
to have a maximal contribution to the total CGB energy flux 
at energies ~1 GeV, this result was a consequence of the "GeV 
excess" in the EGRET measurement of the Milky- Way spectrum 
which recent Fermi observations in mid-Galactic-latitudes have 
not reproduced f Johannesson & Fermi LAT Collaboration 2009; 
lAbdo et al1l2009l ) implying it was likely an instrumental effect. 
In this work, we use a Milky- Way spectrum compatible with no 
GeV excess and find the normal galaxy peak to reside at lower 
energies (see §|2]l. 



tainties and its sensitivity to input parameters and assump- 
tions. 

In this paper, we seek to calculate the expected angular 
correlation of the CGB signature due to 7-ray emitting nor- 
mal star-forming galaxies. We consider two quantities: the 
angular auto-power spectrum of the CGB [C]^) and the 
angular cross-power spectrum between the CGB map and 
some galaxy catalog {CJ^). An advantage of the auto-power 
spectrum analysis is that it can be performed immediately 
after Fermi-LAT has obtained a sufficiently deep all-sky 7- 
ray map, with 7-ray data alone. As such, it does not sufi'er 
from uncertainties introduced through the use of galaxy cat- 
alogs, such as issues of completeness and dust corrections. 
On the other hand, even given the additional uncertainties 
mentioned above, taking the cross-correlation between the 
CGB map and a galaxy catalog provides, as it turns out, a 
better way to detect the normal galaxy angular signature, 
because of the large statistics of the large-scale galaxy sur- 
veys. 

This paper is structured as follows. In § [21 we dis- 
cuss the model we adopt to calculate the contribution of 
normal star-forming galaxies to the 7-ray background. In 
§ [31 we calculate the predicted angular auto-power spec- 
trum form star-forming galaxies, and in § [4] we discuss the 
cross-correlation between the normal-galaxy component of 
the CGB and galaxy catalogs. We summarize and discuss 
our conclusions in § [5l 



2 GAMMA-RAY BACKGROUND FROM 
NORMAL GALAXIES 

We follow PF02 to derive a formulation for the mean CGB 
intensity from normal galaxies. We adopt their assumptions, 
and update our calculation with more recent determinations 
of the cosmic star-formation history and of the Milky- Way 
7-ray spectrum. 

The CGB intensity for photons with energy E (in units 
of photon number per unit area, time, solid angle, and en- 
ergy range) is given by 



I{E) 



47r 



Z 



[{l + z)E,z] 



H(z) 



(1) 



where n.y,com is the comoving 7-ray emissivity density, and 
H{z) is the Hubble function. We assume that the differential 
7-ray luminosity (photons per time per unit energy range) 
simply scales as star-formation rate ^'(z) ^nd gas-mass frac- 
tion ^i{z): 



L-f(E,z) = L~,,MW[E). 

tpMW Mmw 



(2) 



where the quantities with the subscript "MW" represents 
those for the Milky Way. With the comoving number density 
of galaxies iigai, the emissivity is then 



i(i5, z) = L^rigai = L 



7,MW 



(E) 



■^MW m(0) ' 



(3) 



where p*(z) = i/'(2:)ngai is the global star-formation-rate 
density|3 Now, assuming that the sum of gas mass and star 



^ Here we have explicitly assumed, unlike PF02, that most of 
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mass is constant in the typical galaxy, the gas-mass fraction 
is simply given by 



lLdz{dt/dz)p4z) 



(4) 



Thus, given the cosmic history of the star-formation-rate 
density, one could compute ij,{z) by backwards de-evolving 
the present-day Milky- Way gas mass fraction, /xmw- The 
assumption of a total baryonic mass of galaxies staying con- 
stant in time is not necessarily realist ic, as star formation 
is pa rtly fueled by newly accreted gas (jProdanovic fc Fields! 
I2OO8I ). However, the overall effect of the details of the gas 
fraction evolution is relatively small (^factor of two, PF02). 
A more realistic modeling of the evolving gas fraction, in- 
cluding the effects of infall, will be addressed in an upcoming 
publication. 

For the presen t stud y, we adopt a model given by 
iHopk ins fc BeacomI (j2006l ) for the global star-formation- 
rate density as a function of redshift, p*{z). For the 
Milky- Way parameters, following PF02 and references 
therein, we use Vmw ~ 3.2Mq yr~^ and hmw ~ 0.14. 
Lastly, we parametrize the Milky- Way 7-ray luminosity as 
Lj,MwiE) = 1.36X 10=*''(S/600 MeV)""" s"^ MeV"S where 
K = 1.5 for E ^ 600 MeV and k = 2.7 for E > 600 MeV. 
This parametrization comes fro m a broken power- law fit to 
the "GALPROP conventional" (jStrong et al.ll2004h model of 
the energy spectrum of the difi'use Milky- Way 7-ray emission 
(which is compatible with no GeV excess) We set the nor- 
malization by requiring that the energy integral of L^,mw 
above 100 MeV is 2.85 x 10*^ photons s"^ (see PF02 and 
references therein). 

In Fig. [TJ we plot, with the solid line, the 7-ray in- 
tensity E^I{E) from normal galaxies, compared with the 
Sreekumar et al. (1998) determination of the CGB from 
EGRET data. The galaxy contribution appears to be im- 
portant in particular for energies between 50 MeV and 1 
GeV. We point out that due to the shift of the spectral 
break in the Milky- Way diffuse emission spectrum from 850 
MeV (which was the location of the break in EGRET data 
which sufi'ered from the GeV excess) to 600 MeV (the lo- 
cation of the break in GALPROP conventional), the peak 
of the normal galaxy contribution correspondingly shifted 
from ~500 MeV in PF02 to ~250 MeV in Fig. |T] in this 
work. Additionally, the contribution of normal galaxies to 
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Figure 1. The CGB intensity from normal galaxies, compared 
with the Sreekumar et al. (1998) determination of the CGB from 
EGRET data and preliminary Fermi data. 



the CGB declines with energy above 1 GeV faster than it 
did in PF02, as the high-energy slope of the Milky Way spec- 
trum adopted here (2.7) is steeper than the value implied 
by EGRET data (2.4) and adopted by PF02. It is worth 
noting that a preliminary analysis of Fermi data indicates 
that the slope of the CGB spectrum at high energies may 
be substantially steeper (consistent with ~ E~^'*^, see M. 
Ackermann for the LAT ColaboratiorQ) than the EGRET 
measurement. To illustrate this point, in Fig.[T] we also plot 
the preliminary Fermi CGB results. 

In Fig. [21 we plot, with the solid line, the integrand of 
Eq. IT} in units of the integral as a function of redshift at 
E = 300 MeV; this quantity represents the contribution to 
the mean CGB intensity at a given energy from galaxies in a 
specific redshift range. Following the evolution of the cosmic 
star-formation rate, it peaks at z ~ 1 and declines for higher 
redshifts. 



the normal-galaxy 7-ray emissivity at a certain redshift comes 
from galaxies with similar 7-ray properties; it is the properties 
of that typical galaxy that evolve with cosmic time according to 
the product of the cosmic star-formation rate and the gas-mass 
fraction histories. 

^ Note that this is not the latest version of GALPROP that is 
used in the Fermi LAT data analysis and in the comparison with 
Fermi LAT data on the diffuse emission from the Galaxy. How- 
ever, the anisotropy models we are discussing here are not very 
sensitive to the details of the input single-galaxy intensity spec- 
trum. The aspects of the cumulative normal galaxy intensity spec- 
trum that affect our analysis the most are the energy of the peak, 
and the energies above which we are in the power-law tail of the 
spectrum; both these issues can be adequately treated using the 
simple models we adopt here, and for this reason we do not engage 
in a more detailed analysis of the cumulative intensity spectrum. 
We will return to the latter in an upcoming publication. 



3 ANGULAR AUTO-POWER SPECTRUM 

FOR THE CGB FROM NORMAL GALAXIES 

The angular auto-power spectrum of the CGB map due to 
normal galaxies is given by 

where r is the comoving distance and Pgai(fc, z) is the galaxy 
power s pectrum at comovi ng wave number k and redshift 
z (e.g., lAndo et al.l l2007bl ). The multipole £ corresponds 
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Figure 2. Contribution from unit redshift range to the mean 
CGB intensity at 300 MeV (soUd) and angular auto-power spec- 
trum at ^ = 10, 20, 50 and 100 (dotted). 



roughly to the angular scale of 9 = 180° /£. Note that we de- 
fined CJ^ as the variance of the fluctuation from the mean 
intensity in units of steradian. 

The galaxy power spectrum is a well measur ed quantity 
according to the rnodern galaxy surveys (e.g. , I Cole et ahl 
I2OO5I : iMaller et all l2005l : jPercival et all |2007| ). It traces 
the underlying matter power spectrum. To compu t e the 
latter, we adopt the halo-model approach (|Seliakl I2OO0I : 
ICoorav fc ShethI l2002f) with the linear transfer function 
given bv lEisenstein fc Hul l|l999l ). which gives a reasonable 
fit to the galax;y power spectrum with a moderate correc- 
tion for the bias, e.g., bgai = 1.11 (|Afshordi. Loh. fc Straus3 
I2OO4I ). 

In Fig. [31 a), we show the angular auto-power spectrum 
£{£+ \)Cf for E = 300 MeV, as a function of multipole 
£. In the muhipole range between 1 and 10^, £(^-|- l)C7^/27r 
ranges from 10"'' to 10"'^, which is much smaller than the 
case of other sources. For instance, in the case of blazars, 
i(l -\- V)C'P 12-K would be no smaller than ^10"* even at 
large angular scales. In the case of dark-matter annihila- 
tion in the extragalactic halos, it could be as large as 0.1 
at ^ = 10^ (|Ando fc Komatsulbood : lAndo et al. I l2007bl ). or 
even larg er in the case of annihilation in the Milky- Way 
subhalos (|Siegal-Gaskinsll2008l : I Andol 120091 ). In addition to 
the amplitude, the shape of the power spectrum might also 
serve as a diagnostic as it is also different for different source 
populations. 

In Fig. [2j we also show the contribution to the angu- 
lar auto-power spectrum from a given redshift range, i.e.. 



d In CJ^ / dz (the integrand of Eq. (O as a function of redshift 
in units of the integral) for multipoles I = 10, 20, 50, and 
100, which, as we show below, are the observationally rele- 
vant scales. For large angular scales, e.g., £ = 10 and 20, the 
dominant contribution comes from low redshifts mainly be- 
cause of the dependence of the integrand in Eq. ((S)). For 



Figure 3. (a) The angular auto-power (77; solid) and cross-power 
(73; dashed) spectra for E = 300 MeV. The cross-correlation is 
taken with a 2MASS-like galaxy catalog, (b) The relative errors 
of the power spectra after 5-yr all-sky measurement with Fermi- 
LAT. The bin width, A£ = 0.5£, is shown as the arrowed line. 



smaller angular scales, on the other hand, the distribution 
develops a second peak at 2: = 1, reflecting the dependence 
on the cosmic star-formation rate. Therefore, in principle, 
we can probe different redshift ranges by observing C]~' for 
various £, even though these quantities are obtained after the 
redshift information is integrated out. Note that the relative 
contributions of different redshift ranges to C^'' are also dif- 
ferent from those to the mean intensity I{E). In Fig. |4l we 
show the same redshift distribution for C/^ , but focusing on 
the lower-redshift range z < 0.1. Towards larger distances, 
correlations of galaxies are averaged out quickly, and this 
effect is more prominent for large angular-scale modes as 
expected. 

We now examine the auto-correlation detectability with 
the Fermi-ljAT. The la errors for CJ^ measurements are 



given by 



(2£ + l)A£/,ky 



C7^ + 



Cp -l- Cjv 



(6) 



where /sky = f2sky/47r is the fraction of the sky mea- 
sured, ^£ is the bin width for which we use 0.5^, and 
Wi — exp(— ^^af /2) is the window function with the an- 
gular resolution ai, ~ 1?2 for 300 MeV photon. The first 
term represents the cosmic variance, and the second the 
shot noise due to finite statistics of galaxy (Cp) and pho- 
ton (Cjv) counts. The Poisson noise due to galaxies Cp is 
obtained by Eq. ([5]) with replacement Pgai(fc, z) '^gJii 
for the comoving number density of galaxies ngai, we use 
10"^ Mpc"^ and assume it is independent of redshifts. We 
thus obtain Cp = 3.4 x 10~* sr. The Poisson noise due 
to finite photon count is given by Cjv = flsky/N^, where 
A'^^. is the number of photons received from Qsky We esti- 
mate N-y = EI{E)AcaTcfiflsky ~ 5.0 x 10*' /sky, where we 
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Figure 4. The same as Fig. [2] focused on low-redshift range. 



used E = 300 MeV, Ah = 6000 cm^ Teff = ™fov/47r, 
Slfov = 2.4 sr (LAT field of view), and assumed a 5-yr all 
sky survey (T = 5 yr). Thus, we obtain Cjv ~ 2.5 x 10~^ sr, 
which dominates the noise term due to finite galaxy counts 
Cp. Note that the uncertainties calculated here are conser- 
vative, as we have not included the effects of smearing with 
energy within the assumed energy bin (AS ~ i5, with the 
bin extending from E to 2E'). Although the Ct are indepen- 
dent of energy as long as we are in the power-law tail of 
the intensity spectrum, the angular resolution ct;,, and thus 
the uncertainties &Ci, do depend on energy. However, since 
in the spectrum high-energy tail the intensity is decreasing 
with energy as ~ the photons in the energy bin will 

be dominated by the low-energy photons, and the effect of 
energy smearing will be small. In addition, as Oh decreases 
with increasing energy, the inclusion of higher-energy pho- 
tons will result in a decrease of the overall uncertainty. 

In Fig.EIb) we plot, with the solid line, iCj' jCj'' as- 
suming all-sky coverage (/sky = !)• There appears to be a 
sweet spot between ^ sa 5 and 70, where one can claim pos- 
itive detection of galaxy clustering in the CGB with 5-yr 
Fermi data. Below this region, as we have only 2£ + 1 modes 
for fixed ^, C'P cannot be constrained very well (cosmic 
variance). For £'s larger than 100, corresponding to S < 1°, 
the errors become exponentially large because of the limited 
angular resolution of Fermi-'LKJ: . 



4 CROSS-CORRELATION WITH GALAXY 
CATALOG 

We now discuss the cross-correlation between the CGB map 
and the existing galaxy catalogs. The angular cross-power 
spectrum is given by 



CI 



'iTTl{E)Ng 



dz- 



,{[l + z\E,z) dNg 
r^ dz 



Figure 5. Redshift distribution of a 2MASS-like galaxy catalog 
(solid) and angular cross-power spectrum at £ = 10, 20, 50, and 
100 (dotted). The redshift dependence of 10 X the mean CGB 
intensity is shown for comparison (dashed). 



gal 



(7) 



where we define dNg/dz as the redshift distribution of galax- 
ies and Ng is the total number of galaxies of the catalog: 



, dNg 

dz—^. 

dz 



(8) 



The galaxy auto-power spectrum can also be computed with 
these quantities as 



cNI 



dz 



H{z) (dNg 
dz 



l_ 



(9) 



For the present study, we consider a galaxy catalog 
similar to the Two M icron All Sky Surv ey (2MASS) Ex- 
tended Source Catalog ( |jarrett et al.|[200ol ). This is a full sky 
(/sky ~ 1), near infrared survey of galaxies whose median 
redshift is around z ~ 0.1 and total number is Ng ~ 10^. 
The redshift distribution d\nNg/dz is shown as a solid 
curve in Fig. [5 l for w hich we used fitting formula given in 
lAfshordi et airi|2004h . 

Using this galaxy catalog and the CGB emissivity at 
E = 300 MeV, we compute the angular cross-power spec- 
trum £{£ + 1)C]^ /2tt, showing it as a dashed curve in 
Fig. [3ja) . The amplitude of the cross-power is larger than 
that of auto-power by about an order of magnitude, which 
would make the former easier to be detected. In Fig. [S] 
we show contributions from unit redshift ranges to for 
£ = 10, 20, 50, and 100 as dotted curves. Unlike the case 
of auto-power spectrum, the redshift distribution is fairly 
similar for different angular scales, because the galaxy dis- 
tribution dNg/dz has a much sharper peak than dl /dz. 
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The Ig errors of the cross-power spectrum is e stimated 
by re.g.. lZhang fc Beacom|[2004l : ICuoco et al.ll2007l ) 



5C7 



(21 + l)A£/,ky 



Cp + Cj\ 



[Cr + Cm,,) 



1/2 



(10) 



where we use Eq. ([9} for Cf in this expression; Cjv,g = 
^sky/Ng — 1.5 X 10~®/sky sr is the galaxy shot noise. The er- 
rors for CJ^ are plotted as a dashed curve in Fig.|3jb), which 
shows similar prospects to the case of auto-power spectrum, 
for detecting the galaxy clustering in the CGB anisotropy. 
The sweet spot is slightly wider than that for the auto-power 
spectrum. 



5 DISCUSSION AND CONCLUSIONS 

In the calculations for both the mean intensity and 
anisotropy, we assumed that the 7-ray emissivity of all galax- 
ies at the same redshift is the same, rescaled from the emis- 
sivity of the Milky Way using the cosmic star-formation rate 
as well as the gas-mass fraction. This implicitly assumes that 
all the galaxies of interest are Milky- Way-like in their 7-ray 
properties. Although this is clearly not true for all galax- 
ies, what it really amounts to is assuming that most -y-ray 
photons emitted by star-forming galaxies come from Milky- 
Way-like, properly de-evolved sources. This in turn is not an 
unreasonable assumption. Milky- Way-like objects are rich in 
both star formation and gas, so they are expected to be the 
most 7-ray bright among normal star-forming ga laxies of 
the same epoch (see, e.g., iPavlidou fc Fields! I2OOII ). This is 
the reasoning behind taking, as a first approximation, all 
normal galaxies contributing to the CGB to have a single 
luminosity at a given redshift, instead of using a luminosity 
function. 

A notable exception to this general rule is that of 
starburst galaxies, which, depending on the details of cos- 
mic ray confinement, could individually be one to two 
orders of magnitude brighter in 7-rays than the typical 
Milky- Way-like galaxy of the same cosmic epoch, as well 
as have harder energy spectra at high energies (see, e.g., 
[Thompson. Quataert. fc Waxmanll2007l ) . However starburst 
galaxies would be best treated as a distinct source class as 
far as their anisotropy properties are concerned — compared 
to normal star-forming galaxies, the population of star- 
burst galaxies consists of few and bright sources, and the 
anisotropy at small angular scales could be considerably 
stronger, even if the overall contribution of starbursts to 
the CGB is lower. We will return to this issue in a future 
publication. 

Until this point, we have concentrated on photons of 
E = 300 MeV, as the 7-ray energy fiux spectrum due to 
normal galaxies peaks at this energy (see Fig. [l]). We now 
also discuss the results for higher energy photons E = 1 
GeV. At this energy, although the energy flux from star- 
forming galaxies is lower, the LAT effective area increases to 
8000 cm^, and also t he angular resolution improves to ~0?8 
jAtwood et al.ll2009l ). In Fig. [S] we show the same plots as 
Fig.[3]but for E = 1 GeV photons. We find that especially for 
the cross-correlation, the detection prospects are still pretty 
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Figure 6. Same as Fig. |3] for photons of E = 1 GeV. 



good even though the number of photons received would 
decrease. In addition, the multipole range for the detection 
becomes larger, above £ = 100 for C]^ . 

Althoug h we fixed the gala xy bias to be 6gai = 1.11 
throughout (|Afshordi et al.ir2004l ). the relevant value of the 
bias might be different. This is because while 6gai refers to 
all galaxies, we are not interested in elliptical galaxies, faint 
dwarfs, or gas-poor dwarfs, which do not emit significant 
amount of 7-rays. Thus, the bias of star- forming 7-ray- bright 
galaxies might differ from that of all the galaxies as inferred 
from galaxy catalogs. However, our results are not very sen- 
sitive to the value of 6gai as long as the true value is not 
significantly different from 1. 

The emission from individual galaxies we have consid- 
ered is entirely due to the interaction between cosmic rays 
and interstellar gas and light; any contribution from point 
sources within galaxies has not been accounted for. How- 
ever, using 7-ray observations of the Milky Way for guid- 
ance, we expect that the contribution of point sources to 
the total gamma-ray emission of a star-forming galaxy is 
small. First, the relative intensity of the diffuse flux is much 
higher than the total emissio n due to resolved poi nt sources 
in the Milky Way (see, e.g., iHartman et al.lll999i l. Second, 
the good agreement between the Milky Way diffuse emis- 
sion as measured by Fermi LAT and GALPROP, indicates 
that in the Milky Way point sources are not a dominant 
component in the diffuse emission, at least at mid-Galactic 
latitudes. The situation can be very different for early-type 
galaxies with little star formation, in which most 7-ray emis- 
sion would arise from nonthermal processes in older popula- 
tions of stellar remnants, such as millisecond pulsars. How- 
ever the contribution of early-type galaxies to the CGB is 
expected to be very small. 

We also comment on the nonlinear part of the galaxy 
power spectrum. The angular scales where such nonlinear- 



ity becomes important are 



10'' for C7^ and £ 



100 



for C/® . Remembering that the angular resolution of Fermi- 
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LAT for 1 GeV photons corresponds roughly to £ ~ 100, 
the nonhnear part of the power spectrum does not affect 
the relevant result much. The reason why nonlinearity be- 
comes important at lower £ for the cross-power is that with 
the 2MASS-like galaxy catalog the contribution is biased to 
lower redshifts (compare redshift distribution in Figs. [2] and 
[5| that correspond to larger spatial scales for fixed angu- 
lar scales (note k = £/r in the argument of galaxy power 
spectrum -Pgai)- 

In conclusion, motivated by the fact that normal galax- 
ies provide a guaranteed contribution to the CGB and that 
Fermi-LAT has a good sensitivity to measure it, we have 
theoretically computed the CGB angular power spectrum 
due to normal galaxies. We have calculated both the auto- 
power {C]~') and the cross-power {C']^) spectra, using the 
well measured galaxy power spectrum; for the cross-power, 
we correlated the CGB with a 2MASS-like galaxy catalog. 
We found that the amplitude of C]^ is smaller than that 
for other sources such as blazars and dark-matter annihila- 
tion. Still, Fermi-LAT can measure the significant feature of 
the galaxy clustering for the multipole range 10 < ^ < 100 
in about 5 years. The amplitude of the cross-power spec- 
trum C]^ is larger, and the detection prospects are better 
for higher-energy photons. We also found that the redshift 
ranges that contribute to the power spectrum the most are 
different from the case of mean intensity. This feature might 
be helpful in probing the 7-ray luminosity density from nor- 
mal galaxies at various redshift ranges. 
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